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ABSTRACT 



Aims. We investigate the contribution to the formation of type la supernovae of the single (a white dwarf accreting from a non- 
degenerate companion) and double (two merging white dwarfs) degenerate scenario, as well as various aspects of the binary evolution 
process leading to such a progenitor system. 

Methods. We use the combination of a population synthesis code with detailed binary evolution and a galactic chemical evolution 
model to predict the metallicity distribution of G-type dwarfs in the solar neighborhood. Because of the very long lifetime of these 
stars, this distribution is a good indicator of the entire chemical history of a region. By comparing the observed distribution with 
those predicted by assuming different type la supernova progenitors and evolutionary parameters (e.g. concerning mass and angular 
momentum loss and common envelope evolution), it is possible to constrain the possible combinations of assumptions. 
Results. We find that in order to reproduce the observed G-dwarf metallicity distribution, it is absolutely necessary to include both 
the single and double degenerate scenario. The best match is obtained when all merging C-0 white dwarfs contribute to the latter. The 
correspondence is also critically dependent on the assumptions about galaxy and star formation, e.g. the use of the two-infall model 
vs. a constant star formation rate. However, this does not affect the previous conclusion, which is consistent with the results obtained 
by investigating type la supernova delay time distributions in starburst galaxies. 
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1. Introduction 

Type la supernovae (SNe la), the most powerful events produced 
by nature, are not only of critical importance as cosmological 
distance indicators (cf. the 201 1 Nobel Prize in Physics), but also 
for the chemical enrichment of galaxies. A large fraction of all 
iron in the solar neighborhood is known to have been produced 
by these events. It is commonly agreed upon that SNe la origi- 
nate from white dwarfs (WDs) that exceed a critical mass and as 
a result undergo a thermonuclear disruption. Many things remain 
unclear however: whether this happens at the Chandrasekhar 
mass or above or below, whether only carbon-oxygen (C-O) 
WDs qualify, whether the disruption is a deflagration or deto- 
nation, etc. Perhaps most striking ly, even the type of prog eni- 
tor system remains uncertain (see iMaoz & Mannuccill2012L for 
a review). Since the WD needs to accrete in order to reach the 
critical mass, it is obvious that only interacting (i.e. binary or 
multiple) stars can produce such events. However, many differ- 
ent scenarios for the accretion process have been proposed. The 
most popular ones are known as the single degenerate (SD) and 
double degenerate (DD) scenario. In the SD scenario (see e.g. 
[Whelan & Iben 1973; Nomoto 1982), one WD accretes from a 
late main sequence (MS) or red giant (RG) companion. The mass 
transfer rate towards the WD is of critical importance for the sce- 
nario to work: if it is too high, the WD will suffer from a RG-like 
expansion, if it is too low, the WD will burn the aquired mass in 
classical nova outbursts and will never attain the required critical 
mass to explode. Hence, a commonly invoked mechanism to ex- 
plain such a moderate rate is a stabilizing wind from the WD. In 



the DD scenario (see e.g. llben & Tutukovl 1 984 IWebbinkI 1 983) . 
the explosion is the result of the merger of two WDs, com- 
ing together due to the emission of gravitational wave radiation 
(GWR). In the case of unequal masses, the least massive of the 
two will be tidally disrupted and accreted by the more massive 
one, which is thus allowed to reach the critical mass. It has been 
previously suggested that such an event would lead to an accre- 
tion induced collapse resulting in a neutron star (NS), and not a 
SN la. However, Piersanti et al. (2003) showed that the inclusion 
of the effects caused by rotation may well solve this problem. 
Recent studies also consider the possibility of subluminous SNe 
la from equal mass C-O WD mergers (Pakmor et al. 201C t), per- 
haps even below the Chandrasekhar mass (van Kerkwiik et alj 
2010) . In the DD scenario, one thus needs a double WD binary 
with a small orbital separation, so that GWR will sufficiently 
reduce the distance between the stars within the Hubble time. 
While it is conceivable that due to the nature of most scenarios 
the explosion mass is always similar, it is very important to ad- 
dress the question of the explosion mechanism, not in the least 
given the mentioned importance as "standard candles" in cos- 
mology. 

Many methods have been attempted to discriminate be- 
tween both scenarios and hence to determin e which (or both) of 
them generates SNe la in nature (again, see aoz & Mannuccil 
120121 for an overview). It then concerns the study of 
whole populations and of individual explosions, both his- 
toric (e.g. Schaefer & Pagnotta 2012) and contemporary (e.g. 



rChomiuket^l2012l for SN 201 Ife). A promising method on 
a global population scale is the study of the delay time distri 



1 



N. Mennekens et al.: Progenitors of type la supemovae and the metallicity distribution of G-type dwarfs 



bution (DTD), the response function giving the number of SNe 
la as a function of time after an instantaneous starburst. DTDs 
can be observationally measured in passively evolving (ellipti- 
cal) galaxies, and theoretically computed with population syn- 
thesis codes, assuming different SN la progenitor channels. A 
comparison of the two then allows to constrain the formation 
models. 

This paper does not address the physical (im)possibility of 
either formation scenario, but assumes that both of them may 
work as described by their various authors. From this premise, 
it is then investigated what scenario (or both) best theoretically 
reproduces population properties that can be observationally 
tested, also depending on various other evolutionary assump- 
tions and parameters, both stellar and galactic. 

Comparing the observed DTDs of SNe la in passively evolv- 
ing starburst galaxies to those obtained iinder various assump- 
tions with a population synthesis code. iMennekens et alJ (1201 Oh 
(hereafter MIO) showed that the SD scenario alone can not ac- 
count for the amount of SNe la required at more than a few Gyr 
after starburst. Furthermore, it was found that most DD SN la 
progenitors undergo a phase of stable, quasi-conservative Roche 
lobe overflow (RLOF) followed by a common envelope (CE) 
phase, as opposed to two successive CE phases. Finally, stellar 
rotation was proposed to mitigate the discrepancy between the 
absolute number of events predicted and those observed. 

MIO also made a comparison with DTDs previously ob- 
tained by other groups using population systhesis methods. The 
preliminary conclusion was that on many points, a good agree- 
ment was found with most of them, however only under the 
obvious but challenging condition that the assumptions of the 
other groups were reproduced as much as possible. It was al- 
ready hinted that, apart from initial distributions, the most im- 
portant differences in assumptions would be those about mass 
and angular momentum loss, as well as the treatment of CE 
phases. Since then, it has become clear to the population synthe- 
sis community that it is indeed necessary to investigate the differ- 
ences in the predictions made by different groups, and to look for 
the causes thereof. Currently, a joint explanatory paper by four 
groups working on the delay times of SNe la by mea ns of p op- 
ulation synthesis codes is in preparation (Toon en et alJl2013j). I n 
that study, the grou p s behind the s t udies o f iClaevs et al. ( 2010h . 
iRuiter et alj (1201 lb . iToonen et akl (1201 2|) . as well as this one, 
homogenized their codes as much as technically possible, and 
showed that this does indeed cause the predictions to converge 
toward one another. This work has indeed shown that the most 
flagrant discrepancies (which are evident when one compares 
theoretical DTDs, especially SD, produced by different groups 
with their own assumptions) are being caused by differences in 
the assumptions on these three most critical topics (mass loss, 
angular momentum loss, and CE evolution). However, it has also 
become clear that some (more minor) differences cannot be ex- 
plained by different parameterizations, but are the result of in- 
herent differences in the code, i.e. differences in the "single star 
tracks", the evolution tracks of individual binary components. 
This is of course to be expected, given the very different ways 
population codes can b e conceived, e.g. with the analytical for- 
malisms bv iHurlev et al. (2002) vs. the use of full binary evo- 
lution, including the effects of accretion induced full mixing in 
stars that accrete from their companion through an accretion disk 
(see further). The influence of these differences on the macro- 
scopic predictions relevant for this study will be discussed in 
Sect. 3.3. 

This paper intends to confirm and extend the results of MIO 
by means of a method independent of the previous, i.e. not only 



the number of SNe la at a certain time (the DTD), but their last- 
ing legacy on the chemical history of a region. In addition, it 
does not limit itself to the study of passive elliptical galaxies, 
but aims to reproduce the chemical history and SN la rate of the 
solar neighborhood, in the actively evolving Milky Way Galaxy. 
Section 2 summarizes the most important elements contained in 
the used evolution codes, as well as the assumptions about bi- 
nary evolution, SN la progenitors and the chemical evolution 
model. Section 3 then discusses the obtained results, while Sect. 
4 summarizes the main conclusions. 



2. Assumptions 

2. 1. The Brussels population synthesis code 

The Brussels population synthesis code uses as input thousands 
of binary evolu tionary calculations performed in detail with the 
iPaczvnskil (Il967i) based Brussels binary evolution code. The lat- 
ter, under d evelopment for over thre e decades, is extensively 
described in Vanbeverenet al. (1998a). The population code it- 
self is elaborated in De Donder & Vanbeveren (2004). We repeat 
here only those elements that are of particular importance to this 
study. One of the most important consequenses of using detailed 
binary evolution results in the population code is that the effects 
of accretion on the further evolution of the gainer star are treated 
in detail. To do this, it is assumed that when a star accretes 
by direct impact, th i s proc ess occurs following the "snowfall" 
model by iNeo et al" l (fT977l) . When, however, the oribtal separa- 
tion is sufficiently wide to lead to the occurre nce of an accre- 
tion disk, it has been shown by iPacketl (|1981) that the gainer 
star will relatively soon be spun up to (near-)critical rotation. 
Therefore, it is assumed that in such case the gainer star will be 
fully mixed, an event k nown as "ac cretion induced full mixing" 
(see Vanbeveren & De Loore|[T994l) . Especially in this case, the 
further evolution of this star will be markedly different than if 
it would have been unaffected. Therefore, it is certainly not suf- 
ficiently accurate to estimate that the time needed to obtain an 
accreting WD or a double WD is simply equal to the nuclear 
lifetime of the least massive star. The binary evolutionary cal- 
culation s themselves have been published by IVanbeveren et al.l 
(Il998bh . 

It is common in population synthesis codes to model many 
aspects of binary evolution by the use of parameters. This is nec- 
essary in order to be able to follow a large population of sys- 
tems simultaneously, as is required in this study. Apart from ini- 
tial distributions, the assumptions that most affect the results are 
those about mass loss, angular momentum loss, and the treat- 
ment of CE phases. Given the already mentioned wide variety 
of assumptions encountered in the literature, a small elaboration 
and explanation of our assumptions is in order 

The population code starts from a 10^ M© starburst with a 
given binary frequency (the fraction of primaries that have a stel- 
lar companio n) and metallicity Z . Primary star masses Mi are 
drawn from a iKroupa et al.l (Il993i) initial mass function (IMF) 
normalized between 0.1 and 120 Mq. Three different mass ra- 
tio distributions are considered: a flat one (as standard), one fa- 
voring high mass ratios (Garmanv et al. 1980) and one favoring 
small mass ratios (Hogeveen 1992), all normalized to allow sec- 
ondaries wit h a mass M 2 between 0.1 M© and Mi. Reflecting 
the results of lAbtl (Il983l) . the initial orbital period distribution is 
taken to be logarithmically flat, and is normalized between 1 day 
and 10 years. 

Whenever the initially most massive star (subscript 1) in a bi- 
nary fills its Roche lobe, this will initiate a mass transfer phase. 
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Depending on whether this happens during core H burning (case 
A), shell H burning (case B) or after core He burning (case C), 
the nature of this phase and the consequences for the binary sys- 
tem will be different. 

In the event of case A or early case B (known as case Br, 
i.e. at a time when the donor's outer layers are still radiative), 
mass loss will lead to the star retreating within its Roche lobe, 
and therefore be self-mitigating. The result will be a dynami- 
cally stable stream of matter from donor to gainer How much of 
this matter is actually accreted by the gainer star (subscript 2) is 
determined by the quantity fi, defined by 



M2 



M, 



(1) 



In the case of stable RLOF, p is equal to which is a 
common parameter in most population synthesis codes. In our 
standard model, I3„,ax - 1 , signifying that the donor accepts all 
matter lost by the donor We obviously discuss the influence on 
the results of smaller values, and this discussion will indeed turn 
out to support a value close to 1 in reality. 

Apart from changing the masses of both donor and gainer 
(except when fi - 0), another important consequence of mass 
transfer is a change in orbital period and thus separation. In the 
case of conservative mass transfer (fS - I), this change is simply 
given by 



Pi 



,MifM2f 



(2) 



with subscripts i and f standing for initial and final with regard 
to the RLOF-phase. However, if yS < 1 an assumption is required 
about how much angular momentum is carried from the system 
by a given amount of mass loss. This comes down to specifying 
how and where the matter leaves the system. A commonly made 
assumption in population synthesis codes is that mass is lost 
with the specific orbital angular momentum of the gainer star. 
This is however only justified if the matter leaves the system in a 
way that is symmetric with respect to the equatorial plane of the 
gainer (e.g. through an enhanced stellar wind or bipolar jets). As 
there is limited evidence for significant mass loss in such a way 
in intermediate mass non-degenerate stars, our standard model 
is that matter will instead form a non-corotating circumbinary 
ring after passing through the second Lagrangian point L2- In 
this case, the lost angular momentum is the specific orbital an- 
gular momentum of the disk, which is determined by its radius, 
and considerably larger than that of the gainer (since it is fur- 
ther from the center of mass). Hence, in our standard model, a 
given amount of mass loss will lead to a much larger angular 
momentum loss, and thus more systems that merge, than with 
the specific gainer orbital angular momentum loss. Obviously, 
we will also model the latter assumption. This can be done with 
the prescription proposed by Rappaport et al. ( 1983,) (with A the 
orbital separation and P the orbital period): 



y = V^Mi (1 -/3) 



2nA^ 



(3) 



where 77 is a parameter describing the angular momentum loss J. 
Whenever it can be assumed that rj is constant (implying that the 
location and manner through which matter is lost are constant) 
during the considered timeframe, the period evolution is given 
by (see e.g. lPodsiadlowski et al.lll992h : 



Pi 



Pi 
Pi 



forO < j6 < 1: 

/ M[f + M2f 
~\Mii+ M2i 

foryS = 0: 
M[f + M2f 



If 

Mii 



M2i 



Mii + M2i l\M 



(4) 



(5) 



In the case of mas s loss through L2, it was also shown by 
ISoberman et al.l (Il997h that the value of rj is fairly constant, and 
can be well approximated by 77 = 2.3. This value is based on the 
reasoning that smaller circumbinary ring radii (corresponding to 
smaller 77) will lead to fragmentation and fall-back of the ring due 
to tidal forces by an uneven potential. If specific gainer orbital 
angular momentum is assumed, it is straightforward to show that 

Tj - ( iJ^+M2 ) ■ ^^^^ means however, that rj continuously changes 
during mass transfer, and that Eqs. ^ and (|5]l are in principle 
not valid if applied in one step from the situation before to after 
RLOF. It has been checked however that the method of averag- 
ing 77 before and after the mass transfer phase (as is done in the 
population synthesis code) results in virtually the same result as 
fully integrating the period evolution. 

If the mass ratio between the two stars is extreme {q - 
M2/M1 < 0.2 ) at the onset of mass transfer, this will result in 
an instability dParwi ^ II879h . The donor star will be unable to 
extract sufficient angular momentum from the orbit to remain 
in synchronized rotation, resulting in the mass transfer episode 
quickly becoming dynamically unstable. Tidal interaction will 
cause the secondary to spiral into the donor's outer layers, a pro- 
cess that is treated identically to the unstable type of mass trans- 
fer described in the next paragraph. In this case, it is assumed 
that /3 - regardless of /3max- To ensure continuity, a linear inter- 
polation is used between q - 0.2 (J3 - 0) and q = 0.4 (J3 = fimax)- 

When mass transfer starts at a time when the outer layers of 
the donor are deeply convective (case Be or case C), mass loss 
will not cause the donor to shrink, but to expand. Therefore, its 
expansion will not be slowed down, and an unstable mass trans- 
fer phase will ensue. After a short time, the envelope of the donor 
will engulf the other star, and both will rotate within one com- 
mon envelope. The viscous friction of both stellar cores within 
this envelope will cause a reduction of the orbital separation. Part 
of the lost orbital energy will be converted into kinetic energy 
which can be used to expel the common envelope. If this conver- 
sion is sufficiently efficient, the envelope may be ejected in time 
to allow a new binary with smaller orbital separation to emerge. 
If not, both stars will merge. Many different formalisms have 
been proposed to model the energy conversio n. The most popu- 
lar one is the o'-formalism bv lWebbinld (11984^. which assumes a 
conservation of energy to balance the angular momentum before 
and after the CE phase: 



/l^Roche 



MifM2i _ MiiM2i 
2Af 



2Ai 



(6) 



In this formula, /?Roche is the Roche radius of the donor, while 
A is determined by the density profile of this star's outer lay- 
ers, and a is the energy conversion efficiency. The latter two 
values are quite uncertain, and are in population synthesis con- 
text often combined into a single parameter aA. Our standard 
model takes aA - \ (as do most population synthesis stud- 
ies), while other (mainly smaller) values are also considered. 
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iDewi & TaurisI (l2000l) determined values for A for different kinds 
of stars, mostly finding values well below unity. The Brussels 
code allows to assign a different value for aA to each system ac- 
cording to these r esults, and a run in this mode will be discussed . 
In recent studies, iDavis etal J (120101) and .Zorotovic et a 
have determined values for a also resulting in aA up to an order 
of magnitude below 1 . 

On the other hand, iNelemans et al] ( l2000l) propose the y- 
scenario, which conversely assumes a conservation of angular 
momentum to arrive at a balance of energy: 



AJ_ _ AMiotd _ Mii - Mil 



(7) 



the parameter y in this equation having the same function of de- 
noting conversion efficiency as a in the other formalism. Based 
on observational constraints, the authors propose a value of 
y = 1.5. This yields the following orbital separation evolution: 



A, 



Mif + M2i 
Mii -I- M2i 



7 



Mn - Mif 



(8) 



In any case and in either formalism, for surviving systems it 
is assumed that a CE phase is of such short duration (evidenced 
by the fact that it has never been observed with certainty) that the 
gainer is incapable of accreting any matter, resulting in jS = 0. 

It is possible for a donor star that has performed case A 
or B (either Br or Be) RLOF to fill its Roche lobe again for 
a second time later in its evolution, as a helium star. This is 
known as case BB mass transfer (see i De Greve & De Loord 
[l976; Delgado & Thomas 1981). Because the mass transfer rates 
involved are typically much lower than during the first RLOF, it 
is assumed that the non-degenerate gainer will accrete all of the 
matter lost during this episode, i.e. yS = 1. 

At a later point in the evolution, after the originally most 
massive star has become a WD, the originally least massive star 
will also fill its Roche lobe, initiating a mass transfer episode in 
the other direction. Because of the fact that the accreting star is a 
WD with a small surface, and because the mass ratio between the 
two stars is often very large, it is assumed that this will always 
lead to an unstable episode, resulting in a spiral-in. Therefore, it 
is modeled as a CE phase. One exception is when the combina- 
tion of companion mass and orbital period identifies the system 
as a SD SN la progenitor (see next subsection). 

2.2. Type la supernova progenitors 

In this study, mainly the two most popular formation channels 
for SNe la are considered, the single degenerate (SD) and dou- 
ble degenerate (DD) scenario. Various others, includ ing the so 
called core degenerate channel ("Kashi& Soke^lIoilft . various 
delayed detonation models (e.g. Hachisu et al.l 2012l) 7etc. have 
tentatively been included in our population synthesis code, but 
none of them has proven to be able to result in a notable number 
of SNe la at the right time (see Sect. 3). Interestingly, the pro- 
genitor parameter spaces of SD and DD events do not overlap, 
implying that it is possible that both formation channels are at 
work simultaneously. 



2.2.1 . Single degenerate scenario 

For the SD scenario progenitors, we do not explicitly calculate 
the mass accretion on the WD ourselves, but we make use of the 
results of the groups that treated this scenario in detail. When in 



a binary star consisting of a non-degenerate star and a W D, the 
former fills its Roche lobe, iHachisu et alJ (119991 120081) identi- 
fied those zones (dependent on the WD mass) in the (companion 
mass, orbital period)-plane for which the mass transfer rate will 
fall within the naiTow zone in which stable accretion up to the 
critical mass (assumed to be the Chandrasekhar mass of 1 .4 M©) 
is possible. There are two such zones: one for (late) MS compan- 
ions, and one for RG companions, the latter requiring a larger 
orbital period. It should be noted that most population synthesis 
studies fail to repro duce a numer ous amount of systems in the 
latter zone (see e.g. lHan & Podsi adlowski 2004). This is suspi- 
cious, since systems of this nature (though not all necessarily 
SN la progenitors) are indeed quite abundant ly observed, and 
known as the so called symbiotic stars (see e.g. lBelczvriski et akl 
i200 0). The size and location in the (M2, /')-plane of the zone 
with WD+MS progenito rs is strongly depende nt on the assump- 
tions about the WD wind. IHachisu et all (|2008|) propose the mass 
stripping effect, in which a strong wind from the WD blows away 
a significant amount of the matter lost by the non-degenerate 
star, thereby allowing some systems to avoid the CE phase they 
would otherwise have undergone. The strength of this wind is 
regulated by the parameter ci, which determines the amount of 
stripped matter: 



Strip 



dM, 



df 



ci- 



wind 



df 



(9) 



The authors proposed a generic value of ci = 3, while it is 
stated that observations allow values anywhere between (turn- 
ing the effect off) and 10. One of the aims of the current paper 
will be to constrain this parameter further. As an illlustration, a 
ci of 3 about triples the area of the WD-hMS progenitor zone 
compard to ci =0. Finally, the efficiency of the SD scenario 
also critically depends on the metallicity, since in low metallicity 
systems, the W D wind will not be suffic iently strong to regulate 
the mass stream. iKobavashi et al.l(ll998l) calculate d the zones for 
Z=0.004, so that combination with the results of IHachisu et al.l 
( 1999, 2008) for Z=0.02 allows to interpolate for any Z. In the 
population code, whenever a WD + non-degenerate system falls 
within such a SD SN la progenitor zone, it is assumed that this 
system will result in a SN la within a time that is negligible com- 
pared to the previous nuclear timescale of the stars. 

2.2.2. Double degenerate scenario 

Concerning the DD scenario, it is assumed that every merger 
of two WDs, together meeting or exceeding the Chandrasekhar 
limit of 1.4 Mq, will result in a SN la at the time that GWR 
has reduced the distance between both to zero. It should be re- 
marked that this spiral-in timescale is by no means negligible 
compared to the time needed to obtain a double WD. Because 
the exact physics of the explosion are as yet unclear and it is 
widely believed that the scenario only works for the merger of 
two C-O WDs, our computations are also performed with this 
additional restriction. Pursuant to recent speculation on the mat- 
ter, we also investigate the influence on the population assuming 
sub-Chandrasekhar C-O WD mergers to result in SNe la. 

It is of importance to repeat that there are two typical ways 
in which a double WD, eventually resulting in a DD SN la, can 
be formed. To obtain a double WD in an interacting binaryQ the 
primary and secondary need to fill their Roche lobes succes- 
sively, resulting in two mass transfer phases. The first of these 



' In a non-interacting binary, the orbital separation of the double WD 
will be much too large for it to merge within a Hubble time. 
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two phases can either be a stable RLOF event (either conserva- 
tive or non-conservative) or result in an unstable CE evolution. 
As explained in Sect. 2.1, the second phase is assumed to re- 
sult always in a CE phase, since the accretor is a WD. Detailed 
descriptions of either evolution channel as well as typical exam- 
ples are given in MIO. The conclusion is that double WDs having 
gone through two successive CE phases produce DD SNe la al- 
ready after a few hundred Myr, while events having undergone 
a RLOF event have a delay time of up to several Gyr, but re- 
quire the RLOF to be (quasi-)conservative to avoid merging of 
the components. 



2.3. The Brussels galactic evolution model 

The galactic code combines different output sets of the popu- 
lation code with a galaxy formation model to self-consistently 
calculate the star formation rate (SFR) and the subsequent chem- 
ical evolution of a galaxy. To provide sufficient room for in- 
terpolation, the chemical enrichment (including SNe la) by a 
starburst is included on the one hand for a population of sin- 
gle stars, on the other for a population of 100% binaries, and 
this for four different metallicities: Z-0 (obtained by extrapola- 
tion), 0.0002 0.00 2 and 0.02. It concerns a fu ll set of yields as 
di scussed bv | De P onder & Van beverenI (l2004l) . based on those 
of IWooslev & Weaveil (11995 ) but in particular, as motivated in 
the former paper, with the iron core-collapse yields divided by 
a factor of two. Combination with the current SFR, binary frac- 
tion and metallicity thus allows to calculate the current chemical 
ab undances and SN la rate a t any time. As was already found 
bv lDe Donder & V anbeveren ('2004), one of the most important 
ways in which binaries significantly alter the chemical history of 
a galaxy compared to the computations with single stars only, is 
by the chemical enrichment (especially in iron) caused by SNe 
la. Therefore, it is critical that the rate of these events be an inte- 
gral part of every chemical evolution model, which is physically 
sound and consistent with other ingredients, and not just "im- 
ported" from observations, let alone a free parameter. 

Concerning the binary frequency, there are two possible as- 
sumptions: either this value is assumed to be constant in space 
and time, or it is treated as a function of metallicity. The lat- 
ter assumption, a binary frequency that linearly increases with 
Z, does not produce any models that satisfactorily reproduce the 
chemical evolution of the solar neighborhood (as will be shown 
in Sect. 3). Therefore, a constant binary frequency is assumed, 
with a value of 70%. This reasonably high value is required in 
order to produce sufficient SNe la (which can obviously only 
occur in binaries) to attain the observed SN la rate and the nec- 
essary SN la-specific chemical (iron-)enrichment. This is con- 
sistent with the results of MIO, where it was found that a very 
high binary frequency is also required (and not even sufficient) 
to explain the absolute number of SNe la observed in elliptical 
(starburst) galaxies. 

Another important ingredient of any galactic evolution code 
is a galaxy formation model. This describes not only how the 
galaxy itself was formed, but more critically also provides a 
way to calculate the SFR as a function of time. The evolution 
of the SFR critically determines the chemical composition, SN 
la rate, etc. An often used galaxy formation model is the two- 
infall model by Chiappini et al. ( 1997). This model assumes that 
the Milky Way Galaxy was formed by two successive gas in- 
fall phases, the first forming the thick disk (including halo and 



bulge), the second the thin disk. The gas infall rate in either phase 
is given by the following relation: 

±ifll^^ = Ae-"^^ + Be-^'-''-'^'^\ (10) 
at 

In principle, this equation obviously also depends on the dis- 
tance to the galactic center, however since we only consider the 
solar neighborhood we are only concerned with applying it there. 
Ti and T2 are the mass accretion timescales for the halo and disk 
phase, which are taken as 2 Gyr and 7 Gyr (at the position of 
the Sun) respectively. ti„„^ is the time of maximum accretion on 
the thin disk, and is equal to 1 Gyr. The values of the scaling 
parameters A and B are determined by the requirement of the 
galaxy formation model to reproduce the cuiTent observational 
constraints. The SFR is assumed to depend on both the gas sur- 
face density o-g(t) and the total mass surface density cr(f). Chiosj 
( 1980) adapted the prescription of .Talbot & Arnett (.1975.) to the 
infall model: 

\ O-(t) j \o-(t„o,,)l 

where v is the star formation efficiency in Gyr"'. Both this pa- 
rameter and k can be different during the first and second infall 
phase. *I'(f) is set to if 0-^(0 is below 7 Mo/pc". The favored 
model takes values of respectively 2 and 1 Gyr"' for the effi- 
ciencies, and k = 1 .5 in both phases. It should be noted that the 
values of the p arameters described he re and used in this paper 
are taken from iMatteucci et al.l (j2009) and differ from those in 
the original two-infall model by Chiappini^LalJ ([1997) which 
were used in previous papers by the Brussels group. Notably, 
the 'new' values result in a SFR that is much less concentrated 
at early times, and instead is more flattened out. 

A second possible galaxy formation model that is considered 
in this paper is simply a constant SFR. The onset of star forma- 
tion is assumed to be delayed by 0.5 Gyr, which is the required 
time to allow the surface gas density to rise to 7 M^/pc^. After 
that, the SFR is constant with its value determined by the cur- 
rently observed gas and star densities. This turns out to be about 
4 Mo/pc^Gyr, resulting in a total star formation roughly equal to 
the case of the two-infall model. 

In both cases, the infalling primordial gas has a composition 
of X=0.76, Y=0.24 and Z^Q. 

Since a similar study bv iDe Donder & Van beveren (2004|), 
studies of the same kind have been performed by Greg gio et al] 
(l2008h . IMatteucci et all (|2009|) (albeit both not with internaUy 
computed SN la DTDs, but with various adopted ones) and 
Kobavashi & Nakasato (201 1) (albeit for the SD scenario only). 
The former three concluded that the best match between predic- 
tion and observation was found when both the SD and DD model 
were combined. 

2.4. Observational galactic parameters 

A critical prediction made by galactic evolution models is the 
abundance of various chemical elements. These are usually ex- 
pressed logarithmically, relative to hydrogen (H), and normal- 
ized to the solar abundance. In this context, the most impor- 
tant abundance is that of the element iron (Fe). For the chem- 
ical yields of a single SN la event, those of the W7 model by 
Iwamoto et al. ( 1999) are assumed, both for SD and DD progen- 
itors. Most importantly, this includes a yield of 0.626 M© Fe per 
SN la. 
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The age-metalHcity relation (AMR) is the evolution of 
[Fe/H] (= log(Fe/H) - log(Fe/H)0) as a function of time. Stars 
that are formed at a certain time are obviously assumed to con- 
tain the metallicity of that instant. Therefore, if these stars are 
still observable today, and their age can be determined, they pro- 
vide a way to reconstruct the chemical history of the Galaxy. 
Stars that are very suited for this purpose are the G-type dwarfs 
(O.8OM0 < M < I.OSMq), since they have an extremely long 
lifetime of up to the age of the Galaxy itself. Therefore, the 
metallicity distribution of these local G-type dwarfs, i.e. the rel- 
ative number of observed G-dwarfs with respective [Fe/H], is in- 
dicative of the entire chemical history of the solar neighborhood. 
Since the DTDs of SD and DD SN la progenitors are markedly 
different, this chemical history should be critically influenced by 
which scenario is at work. 

The observational G-dwarf metallicity distri butions that will 
be use d for comparison in this paper are those bv lHolmberg et alJ 
(l2007h . Two distributions are derived in this paper, one for 
a spherical solar neighborhood and one for a cilindrical one. 
Since in our galactic code the solar neighborhood is defined 
as a cilindric region of 1 kpc around the Sun, mainly the lat- 
ter will be used. It should be noted that a m ore recent analy- 
sis of the Geneva-Copenhagen survey data by iCasa^rande et alJ 
(l20Tll) results in a slightly different G-dwarf metallicity distri- 
bution, which has the same morphological shape but peaks at a 
slightly higher value o f [Fe/H]. As the differ ence (~ 0.1 dex) 
is not too large, and the lHolmberg et all ( l2007l) distribution pro- 
vides the advantage of having been explicitly calculated also for 
a cilindrical solar neighborhood, the latter will be used for com- 
parison. Additionally, if theoretical models are ruled out by be- 
ing unable to produce enough Fe, these conclusions will not be 
affected if an even higher [Fe/H] must be matched. 

While the G-dwarf metallicity distribution seems quite ro- 
bust (those alrea dy mentioned are quite similar to the ones by 
iLee et aL I (l20Tl . and also numerous others found in the liter- 
ature over the past decade), the same is not true for the AMR. 
As will be elaborated on in Sect. 3.2, the recent literature con- 
tains a wide variety of such observed relations, depending on the 
method of their construction. Some of these AMRs show a ris- 
ing trend in [Fe/H] as a function of time, while others do not. 
Also the width of the AMR (i.e. the extent of scatter in [Fe/H] 
at a given age) and its possible origins is a matter of debate. For 
this reason, our comparison between prediction and observation 
will focus mainly on the G-dwarf metallicity distributions them- 
selves, and less on the AMR. 

A chemical evolution model also needs to explain the evolu- 
tion of the abundance of other chemical elements, in particular 
[O/Fe] and [C/Fe] vs. [Fe/H]. Especially [C/Fe] is influenced by 
binary effects other than SNe la. 

2.4.1. Comparison method 

The current age of the Galaxy is taken to be 13.2 Gyr, the time 
after which the theoretically obtained star and gas surface densi- 
ties are com pared to the curren tly observed ones. The latter are 
taken from Calura etal ](|201(ft and references therein, yielding 
37.5 ± 10 Mopc"^ for the steflar surface density and 10.5 + 3.5 
Mopc"^ for the gas surface density. Other parameter-sensitive 
output that can be compared to current values are a SFR of 3.5 
+ 1.5 M(7)PC ~^Gyr~' (according to the same source's analysis of 
iRan a 1991) and a local SN la rate of 0.003 ± 0.002 pc ^Gyi- i 
(according to the analysis of Cappellaro 1996). The latter corre- 
sponds to 3 + 2 SNe la per millennium in the Milky Way Galaxy. 
Only when these values are reproduced to within the observa- 



tional uncertainty is a model considered as possibly valid. The 
final distributions that then need to be checked, the actual goal 
of the study, are the C and O histories and, most importantly, the 
G-dwarf metallicity distribution. 

Using the same iKroupa et al.l (1199 3*) IMF as the population 
code, the galactic code calculates from the SFR how many G- 
dwarfs are being born at any given time. These are obviously 
assumed to contain the metallicity that is prevalent at the time 
of their birth. Taking into account the mass-dependent lifetime 
of each G-dwarf, it is then calculated which fraction of those 
born at each timestep is still observable today. The [Fe/H]- 
distribution of these remaining stars is then plotted in a his- 
togram with bin size 0. 1 dex, ready to be compared to observa- 
tions. Theoretically, every star born at the same time thus has the 
same Fe-content. As mentioned before, however, observations 
indicate that this may not be the cas e in r eality (Ramirez et d] 
l2007tfHolmberg et al.i2007l:ICasagrande et al.l20rih . It may thus 
prove to be necessary to supplement the theoretically obtained 
[Fe/H] at each moment with a random (positive or negative) 
deviation, e.g. distributed following a certain standard devia- 
tion. A possible reason for this "intrinsic scatter" in the AMR 
is the r adial migration of stars over large galactocentric dis- 
tances (i Roskar et alj l2008l) . The possible implications of this 
open question for the results of this study will be discussed in 
Sect. 3.2. 

It should be noted that the early [C/Fe] and [O/Fe] vs. [Fe/H] 
evolution is extremely sensitive to the assumptions about direct 
black hole formation, i.e. whether and how much SN ejecta are 
expelled by massive stars prior to collapse, and what the lower 
mass limit for this process is. The current calculations have been 
performed under the assumption that stars ab ove 40 Mp will 
form a black hole with the hypernova yields of iNakamura et al.] 
(i2001[) . the total ejected mass so chosen that 4 M© of O is 
expelled per event (see De Donder & Vanbeveren 2003). Other 
choices do not more than marginally influence the C and O his- 
tories at [Fe/H] > -1.5, nor the G-dwarf metallicity distribution. 
As a remark, the solar reference values used in the determination 
of [Fe/H], [C/Fe] and [O/Fe] are the actual solar abundances, not 
those predicted by the model. Choosing the latter would obvi- 
ously force the curves through the origin. 



3. Results and discussion 

Anticipating the results in this section, we find that reproduc- 
ing the observed DTD is a necessary but not sufficient condition 
for a certain combination of assumptions to satisfactorily repro- 
duce the observed G-dwarf metallicity distribution. Therefore, 
the first test will again focus on those DTDs. Only when the 
DTD is in agreement with the observed ones can there be hope 
of reproducing the observed G-dwarf metallicity distribtution in 
a second step. A good DTD is however no guarantee for suc- 
cess, not only because of galactic evolutionary considerations, 
but also because of the limited timeframe during which the the- 
oretical DTD can be compared with observation: at early times 
no observational test for the DTD is available, while this part 
obviously also contributes to the galactic chemical evolution. 
Hence there will be a considerable difference between the G- 
dwarf metallicity distributions obtained for the following fic- 
tional DTDs: (a) one matching the observational data points and 
with a late onset just in time to match the earliest such point 
(at 175 Myr); (b) one matching the observational data points but 
with a very early onset of 50 Myr; (c) one matching the obser- 
vational data points and with an intermediate onset of 85 Myr 
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(resulting in the best match for the G-dwarf metallicity distribu- 
tion). 

Additionally, whether a non-fictional model with satisfac- 
tory DTD will also reproduce the observed G-dwarf metallicity 
distribution critically depends on the Z-sensitivity of the model. 
The DD channel is relatively independent of metallicity, its efli- 
ciency even increasing a bit with lower Z. However, the number 
of events produced by the SD channel is decimated when Z is 
lowered from 0.02 to 0.002, and reduces to zero for even lower 
Z. This is of course a result of the WD wind no longer being able 
to stabilize the mass stream (see Sect. 2.2.1). This means that a 
model which solely relies on the SD channel to produce SNe 
la with short delays, will only be able to provide these events 
(and their Fe enrichment) after the metallicity has already been 
enhanced. Hence, the early evolution of the galaxy will be char- 
acterized by Fe-deficiency. 

3. 1. Delay time and G-dwarf metallicity distribution 

In this subsection we will first establish our standard model of 
parameters and other assumptions, as also used in MIO. It will be 
shown that this model, with neither progenitor scenario nor with 
both combined, is satisfactory, even with the suggestions for im- 
provement made in that paper. Thus, we will investigate what 
changes in the assumptions would be required to mitigate the 
discrepancy between prediction and observation. Subsequently, 
the following will be studied: 

- the influence of angular momentum assumptions (a vs. y- 
scenario for CE, values of parameters, modeling of angular 
momentum loss) 

- whether inclusion of any non-traditional SN la progenitor 
scenarios can significantly change the picture 

- the influence of uncertainties intrinsic to the tradiational pro- 
genitor scenarios (what is the best value for ci in SD sys- 
tems, are DD SNe la restricted to C-O WD mergers, do sub- 
Chandrasekhar C-O WD mergers also explode?) 

- all along, also the influence of galactic evolutionary param- 
eters and assumptions (e.g. galaxy formation model, binary 
frequency type and value) will be studied. 

The DTDs obtained with the population code for a sin- 
gle starburst consisting of 100% binaries are shown in Fig. [T] 
Represented are the SD DTD for /3 max = 1 and ci = 3, as well 
as the DD DTDs for /3max - 1, respective l y 0. Su perimposed are 
the observational DTDs of 'Totani et al.' ('2008) (including one 
data point by Man nucci et a l. 2005) and Maoz et al. (2012). The 
SD DTD for = barely differs from the one shown. Other 
assumptions are Z=0.02, aA = I and a flat mass ratio distri- 
bution. As already concluded by MIO, independent of the latter 
choices, the SD DTD drops away too fast and too soon to keep 
matching the observations after a few Gyr, even with values for 
ci up to 10. For /3max - 1 the DD distribution matches the ob- 
servations in morphological shape, but lies a factor 2-3 too low 
in absolute value at the 1 1 Gyr point. The drastic decrease in 
the DD DTD for /3max = means two things: firstly, the events 
which disappear in the latter case (about 80% of all DD SNe 
la) have all gone through a stable RLOF event and thus not two 
CEs; secondly, in order to be able to reproduce the morpholog- 
ical shape of the observations, /Jmov needs to be close to 1, i.e. 
> 0.9. Concerning the SD events, as already touched upon in 
Sect. 2.2.1, we also find very few WD+RG progenitors, most 
events occurring in WD+MS systems. Also shown in the figure 
is the combination of both SD and DD, albeit with ci = 1 and 



assuming a 10% increase in convective core mass due to rotation 
(as proposed by MIO), the latter mainly enhancing the DD chan- 
nel. This combination matches the observed DTDs both in shape 
and number, but as will be seen later this does not automatically 
mean that inserting this DTD in the galactic model will yield a 
satisfactory galactic chemical evolution outcome. 

Fig. |2] shows the G-dwarf metallicity distribution obtained 
with only the SD and DD channel respectively. The SD rates 
have been calculated with ci = 3, a moderately high value which 
results in rates only marginally smaller than if the maximum 
value of ci - 10 would be chosen. While the DD DTD is the 
only one still compatible with the observations in morphological 
shape, it is obvious that the G-dwarf metallicity distribution does 
not correspond to observations. The same is true to an even more 
dramatic extent for the SD channel. In none of the performed 
simulations, with none of the stellar or galactic evolutionary pa- 
rameter choices, we obtain a situation in which either scenario 
alone succeeds in even marginally reproducing the observed G- 
dwarf metallicity distribution. 

As does Fig. [T] MIO showed that when a 10% increase in 
convective core mass is assumed, the combined SD + DD model 
matches the observational DTDs well. However, it is obvious 
from Fig. |2]that this is not true for the corresponding G-dwarf 
metallicity distribution. This shows that the distribution peaks at 
too low [Fe/H] values, indicating that not enough Fe is formed 
early on in the galactic evolution. This conclusion stands inde- 
pendent of the chosen value for ci. The reason is that the early 
part of the DTD, prior to the first observational test, does not 
provide enough Fe pollution soon after the starburst. These early 
events are caused through the SD channel (which is not much af- 
fected by the convective core mass increase and does not work 
for low Z anyway), as well as through WD mergers akeady mas- 
sive enough without the increase. As will be shown later, if early 
events would also increase by a factor of 2-3, the G-dwarf metal- 
licity distribution would match observations. Obviously, the dis- 
tribution for the combined SD + DD model but without convec- 
tive core mass increase (also shown in Fig. 2) results in an even 
poorer match with observations. 

As an academic exercise. Fig. |3] shows the G-dwarf metal- 
licity distribution which is obtained when the standard model 
SD and DD channels (from Fig. [TJ are combined, but multi- 
plied by a factor of 2.5 over the entire time rang^ In that case, 
the DTD perfectly matches the observations. When the standard 
two-infall model is used, this results in a G-dwarf metallicity 
distribution that peaks in the correct location, but with much less 
spread than the observations. On the other hand, when a constant 
SFR is assumed, the predicted G-dwarf metallicity distribution 
shows a good match with the observed one. As a side note, it 
is remarkable that the one obtained with the two-infall model 
shows a better resemblance with the observational distribution 
obtained for a spherical solar neighborhood as is used for com- 
parison by Kobavashi & Nakasato (201 1). 

The DTDs and the resulting G-dwarf metallicity distribution 
were also calculated using the y-scenario for CE evolution in 
the first mass transfer phase. The results are shown in Fig. |4] 
and |5] respectively. The same figures also show the results ob- 
tained with the a-scenario for CE evolution, but with the aA pa- 
rameter not set to the same fixed value for all stars but instead 
ca lcultated separa tely for each system according to the results 
of iDewi & Taurisl (2000), resulting in values for this parameter 



' To allow for straightforward comparison, the same factor of 2.5 will 
be used in several subsequent models that otherwise underestimate the 
observed Fe-abundances. 
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Fig. 1. DTDs obtained with the DD scenario for /3max = 1 (solid) and j6,„(u = (dashed), as well as with the SD scenario for /3max - 1 
and ci = 3 (dotted). Combined SD + DD mo del DTD with /3max - 1 and ci = 1 , but with a conv ective core mass increase of 10% 
(dashed-dotted). Observational data points bv lTotani et alJ (l2008h (circles) and lMaoz et alJ (l2012h (squares). 
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Fig. 2. G-dwarf metallicity distributions obtained with the combined SD (with ci = 3) + DD model either without (solid) or with 
(dotted) 10% convective core mass increase, as well as for the SD (with ci - 3, da shed) respectively DD (dashed-dotted) scenario 
alone. Observational data for a cilindrical solar neighborhood (white histogram) bv lHolmberg et al] (l2007h . 



much smaller than one, on average a A ^ 0.25. It can be sur- 
mized that for both alternative CE treatments the conclusions 
are similar as for the standard ff-scenario: in order to reproduce 
the observed G-dwarf metallicity distribution, the absolute num- 
ber of SN la events needs to be multiplied by 2.5. This then also 
results in a DTD matching the obser ved one. Only when using 
the C/l-values by Dewi & Tauris (2000) does this multiplication 
result in a slight Fe-overabundance, indicating that a lower mul- 
tiplication factor is needed there for the best match. Also notable 
is that when using the y-scenario, the morphological shape of the 
G-dwarf metallicity distribution is acceptable for both the two- 
infall model and the constant SFR. 

Use of the specific gainer orbital angular momentum loss as- 
sumption instead of the L2 one (see Sect. 2.1) results for fi„,ax - 1 
in SD and DD SN la rates changing by not more than some 10%, 
nor is the morphological shape of the DTD severely altered. The 
same is thus also true for the resulting G-dwarf metallicity distri- 
butions and the conclusions from those. One possibly important 



difference is that with this angular momentum loss assumption, 
in the case of y6,„fli - 0, the distributions are not greatly influ- 
enced compared to j = 1 , owing to the very limited angular 
momentum loss. Thus, much lower values of (3 than with the Li 
assumption still yield DTDs compatible in morphological shape 
with the observations. 

When the standard SD + DD model is considered, but with 
DD SNe la restricted to the merger of two C-O WDs (and still 
exceeding the Chandrasekhar limit), the total number of SNe 
la goes down by some 20%, but this does not compromise the 
morphological shape of the DTD. The corresponding G-dwarf 
metallicity distribution, again for the DTD multiplied by a fac- 
tor of 2.5, is also shown in Fig.|5] It demonstrates that to obtain 
a satisfactory distribtution, the multiplication factor then merely 
has to be increased slightly. 

One of the non-standard SN la formation sc enarios that has 
been tested is the core-degenerate (CD) channel dKashi & Soked 
201 li) . For this, it is assumed that the merger of a WD 
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Fig. 3. G-dwarf metallicity distributions obtained with the combined SD (with ci = 1) + DD model, with SN la rates multiplied by 
a factor of 2.5, for a constant SFR (solid) and for the tw o-infall model (dotted). Observational data for a spherical (gray histogram) 
and cilindrical (white histogram) solar neighborhood bv lHolmberg et al.] (l2007h . 
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Fig. 4. DTDs obtained with the y-scenario for CE-evoluti on throug h the DD (solid) and SD (dotted) channel. DTDs obtained with 
the Qf-scenario, using the parameter v alues determined bylDewi & Tauris (2000), th rough the DD (dashed) and SD (dashed-dotted) 
channel. Observational data points bv lTotani et al.l (l2008l) (circles) and Maoz et all (1201 2h (squares). 



with the core of a non-degenerate star, together exceeding the 
Chandrasekhar mass, will result in a SN la. This merger can 
happen at any time between the formation of the first and sec- 
ond WD, but most U kely during a CE phase. According to 
lllkov & Sokerl (l2012ah . the SN la does not necessarily take place 
at the time of merger: if they have a rapid rate of rigid rota- 
tion, the explosion of merger products with a mass below 1 .48 
Mq can be delayed by a significant time (up to 10 Gyr), un- 
til magnetic braking has provided enough downspin. This can 
not only deliver the events with long delay time, but can also 
explain the absence of H in the spectrum of the explosion. It 
should also be noted that the CD channel is relatively indepen- 
dent of the aA parameter of CE evolution. If this product is set 
to O.I instead of 1, the CD rate changes by only -15%, compared 
to -95% and +80% for SD and DD respectively. Also the cor- 
responding change in morphological shape of the distribution 
is much smaller than for those two. However, as was akeady 



shown by MIO (without making reference to SNe la) the num- 
ber of WD + non-degenerate mergers is extremely large imme- 
diately after starburst, but then drops very fast (with the number 
of such events sharply falling by several orders of magnitude af- 
ter 0.3 Gyr). Although the total number of these potential CD 
ev ents matches the total number of SNe la quite well (as noted 
by lllkov & Soker|[2012bh . we thus find that the number of early 
events (< 0.3 Gyr) lies much too high to match the observed 
DTD, while the number of late events (> 0.4 Gyr) lies much too 
low. The mentioned delay between merger and explosion does 
not solve this problem, as it can only apply to systems with a 
total mass < 1.48 M©, which represent only 10% of the total. 
Also the G-dwarf metallicity distribution obtained under the CD 
assumption is unable to match the observed one: owing to the 
extremely high number of events (up to > 10"* Gyr"') at early 
times, it results in a strong Fe-overproduction regardless of other 



9 



N. Mennekens et al.: Progenitors of type la supemovae and the metallicity distribution of G-type dwarfs 




[Fe/H] 

Fig. 5. G-dwarf metallicity distributions obtained with the combin ed SD (with c\ = 3) + DP model, all with SN la rates multiplied 
by a factor of 2.5, for the a-formalism with parameter values by iDewi & TaurisI (l2000h (solid), for the y-formalism and either a 
constant SFR (dotted) or the two-infall model (dashed), as well as for the standard c-formalism but with only super-Chandra CO- 
C O WD mergers resultin g in DD SNe la (dashed-dotted). Observational data for a cilindrical solar neighborhood (white histogram) 
bv lHolmberg et all (|2007[) . 
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Fig. 6. DTD obtained with the assumption that all double C-O WD mergers result in SNe la, through the DD (solid) channel. 
Combined SD (with ci = 1) + DD model DTD under the same ass umption, bu t with a convective core mass increase of 10% 
(dotted). Observational data points bv lTotani et al.l ( l2008l) (circles) and lMaoz et al.l ( 1201 2|) (squares). 



assumptions. These conclusions also stand under the additional 
restriction that the WD involved must be of the C-O type. 

For the same reason, the delayed detonation SD model pro- 
posed by Hachisu et al. (2012) is unable to provide a significant 
increase in SNe la with long delay in our simulations. This sce- 
nario assumes that rotating WDs can accrete matter until highly 
exceeding the Chandrasekhar limit, not exploding until the star 
has spun down. However, only for the small minority of rigidly 
rotating WDs with a mass below 1 .5 Mq can this cause a signif- 
icant delay. 

Recently there has been much speculation about the possi- 
bility of su b-Chandrasekhar C-O W D mergers exploding as SNe 
la (see e.g.'Badenes & Maozl l20T2h . Fig. [f] shows the DD DTD 
obtained under the assumption that all C-O WD mergers, regard- 
less of their total mass, indeed result in such an event. This DTD 
(in combination with the SD channel from Fig.[T]i shows a very 



good agreement with the observational ones. However, as Fig. 
Q shows, the obtained G-dwarf metallicity distribution shows 
not enough Fe enrichment, even though the discrepancy is not 
too large. The reason for this is again that at the very start of 
the DTD, before any comparison with observation is available, 
the SN la rate is low, especially for Z<0.02 (DTDs not shown). 
However, these very early events (< 0.1 Gyr) seem to be re- 
quired in order to cause enough Fe pollution at very low Z (and 
thus [Fe/H]). Of course, this Fe might also have another source 
than SNe la. A simulation was run including the yields of very 
m assive stars (140 Mp < M < 260 Mq) at Z < 10"^ as computed 



bv lHeeer & Woosle vr(l2002h. which liberate u p to 40 Mq of Fe 
each (see also Vanbeveren & De Dondeill2006h . Because the du- 
ration of this low metallicity phase is so short, the Fe produced 
during this era only greatly influences the chemical evolution at 
that time itself, and does not leave a lasting legacy once other 
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Fig. 7. G-dwarf metallicity distributions obtained with the combined SD + DD model, with ci - 3 under the assumption that all 
double C-O WD mergers result in SNe la (solid), with ci - 1 and the same assumption but for a convective core mass increase of 
10% for either a constant SFR (dotted) or the two-infall model (dashed), as well as for a constant SFR but with SN la yields scaled 
to their total explosion m ass if < 1.4 Mq (dashed-dotted). Observational data for a cilindrical solar neighborhood (white histogram) 
bv lHolmberg et all (I2007I) . 
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Fig. 8. G-dwarf metallicity distributions obtained with the combined SD (with ci - I) + DD model, with all double C-O WD 
mergers resulting in SNe la and a convective core mass increase of 10%, for a constant binary frequency of 35% (solid) and for a 
binary freqency increa sing linearly with Z (dotted). Observational data for a cilindrical solar neighborhood (white histogram) by 
iHolmberg et"aLl(l2007h . 





sources of Fe (i.e. mainly the SNe la) turn on and s oon drown it 
out. Th e origin of this early Fe thus remains unclear IWang et aP 
(I2009h propose a He-star channel which they claim leads to very 
early SNe la through the SD channel. When introduced into our 
population code however, this channel seems to produce only 
relatively late SNe la, and not many compared to the other chan- 
nels. This seems logical, since in the SD channel one always 
needs the secondary to fill its Roche lobe towards the WD. As 
the maximum allowed secondary mass for this mass transfer to 
be stable is around 3 Mq, it is obvious that this can take a while. 
One hypothetical solution to the situation is when one does not 
restrict the DD SN la explosions to C-O WDs. If, in addition to 
assuming that all C-O WD mergers explode (regardless of total 
mass), one also assumes that all super-Chandrasekhar mergers 
explode (regardless of composition), then the early-time events 



(starting as early as 20 Myr) are provided by DD mergers where 
the more massive star is an ONeMg WD. Without compromising 
the shape and number of the DTD at times for which observa- 
tions are available, the obtained G-dwarf metallicity distribution 
then shows an excellent match with the observations. 

A much less artificial model, that also reproduces the ob- 
served G-dwarf metallicity distribution, is when the assumption 
about all C-O WD mergers exploding as SNe la is combined 
with the assumption about the 10% convective core mass in- 
crease. The then obtained DTD is also shown in Fig. |6] and 
perfectly matches the observations, without multiplication fac- 
tor The same is true for the corresponding G-dwarf metallicity 
distribution shown in Fig. [7] This figure shows three different 
distributions obtained under these assumptions. Firstly for the 
flat SFR vs. the two-infall model, both for ci - I. This is the 
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time (Gyr) 

Fig. 9. AMR with the combined SD (with ci = 1) + DD model, with all double C-O WD mergers resulting in SNe la and a 
co nvective core mass in crease of 10%, for a constant SFR (black) and for the two-infall model (gray). Observati onal data points 
by 'Rami rez et alj (12007*) (dots) and for the Sun (open cricle), as well as zones (shaded) in which the data points bv lHolmberg et al] 
(2007) (lower) and Casagrande et al. (2011) (upper) lie, but see text for important note on their representativeness. At the lower 
right, the typical observational error bar for center of the plot (increasing towards the left) is shown. 

6 1 




time (Gyr) 

Fig. 10. SN la rate with the combined SD (with ci - 1) + DD model, with all double C-O WD mergers resulting in SNe la and a 
convective core mass increase of 10%, for a constant SFR (black) and for the two-infall model (gray). Currently observed rate is 
indicated by error bar. 



only model in which for a flat SFR stellar surface density, gas 
surface density and current SN la rate are within observational 
constraints at the same time. The model with ci = 3 results in a 
slightly too high SN la rate (the same was true in the case of the 
combined SD + DD model multiplied by a factor of 2.5). The 
third distribution shown is also for a flat SFR, but with all DD 
SN la yields scaled to their total explosion mass (if below 1 .4 
Mo), i.e. taking into account the fact that a sub-Chandrasekhar 
DD SN la may well produce less Fe (and other yields) than a 
Chandrasekhar mass one. The figure shows however that this 
open question does not critically affect the metallicity distribu- 
tion. Figure |2] shows a slight Fe-overproduction when using the 
SFR from the two-infall model. In that case, the assumption that 
all C-O WD mergers result in DD SNe la yields a better G- 
dwarf metallicity distribution without the additional assumption 
of a 10% convective core mass increase than with it. However, it 



then only concerns the location of the peak of the distribution, its 
morphological shape is still reproduced better with the flat SFR. 

Even in the case of the best model, the predicted G-dwarf 
metallicity distribution shows too little (i.e. no) stars with very 
high metallicity, and too many with very low metallicity. This is 
the case for all reasonably simple theoretical galactic chemical 
evolution models, and is long known as the G-dwarf problem. A 
possible way to mitigate the situation is by assuming that stars 
are not necessarily observed at the location of their birth (and 
for whic h their [Fe/HI is thu s representative). Such a migration 
scenario dRoskar et al.ll2008l) allows to explain the observational 
presence of stars with higher metallicity than can be produced 
locally according to the model. Additionally, it allows to repro- 
duce to a certain extent the intrinsic scatter in the AMR that is 
observed (i.e. the fact that not all stars born at the same time have 
the same [Fe/H] in reality, see also next subsection). 
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Fig. 11. [C/Fe] vs. [Fe/H] with the combined SD (with ci = 1) + DD model, with all double C-O WD mergers resulting in SNe 
la and a convective core mass increase of 10%, for a constant SFR ( black) and for the two-infall model (gray). Observational data 
by IChiappini etal] (l2003h . iReddv et all (l2003l) . iBensbv & FeltzingI (|2006|) and iFabbian et all (|2009|) (dots) and for the Sun (open 
cricle). 
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Fig. 12. [O/Fe] vs. [Fe/H] with the combined SD (with ci - 1) + DD model, with all double C-O WD mergers resulting in SNe la 
and a convective core mass increa se of 1 0%, for a constan t SFR (black) and for the two-infall model (gray). Observational data by 
[Gratton et al. (2003) . Reddv et all (l2003h and lCalura etaP (1201 Ol) (dots) and for the Sun (open cricle). 



As an illustration of the previously claimed requirement to 
have a high and constant binary frequency, Fig.[8]shows the same 
best model but for a binary frequency of 35% (half the standard 
value), respectively one varying linearly with Z. It is obvious that 
the G-dwarf metallicity distribution in both cases shows a very 
strong underproduction of Fe. 

3.2. Other galactic evolutionary constraints 

Figure|9]shows the AMR for the model which best approximates 
the observed G-dwarf metallicity distribution, i.e. the SD (with 
ci - 1) + DD model but with all C-O WD mergers exploding as 
SN la and a convective core mass increase of 10%. This is done 
both for a constant SFR and for the SFR following from the two- 
infall model. It should be noted that distributions for other (rea- 
sonably) satisfactory models, such as the same model without 



the core mass increase or the standard SD + DD model with SN 
la rates (artificially) enhanced by a factor of 2.5, look very sim- 
ilar When the pre dicted AMR is comp ared to t he observational 
ones obtained by Holmberg et alj (l2007ft and iCasagrande et al.l 
(20 11), this is not very instructive. The latter namely suggest an 
AMR that is almost flat, but has a very significant and intrinsic 
scatter As a result of the latter, all theoretical models produced 
in this study lie within the observationally populated zone, but 
none of them is able to reproduce an AMR with this general (lack 
of) shape. This is because the infalling pristine gas is unable to 
fully compensate for the larg e Fe-enrichment caus ed by S Ne la. 
If the iHolmberg etakl (l200'7t) or ICasagrande et al.1 (1201 lb AMR 
is indeed representative for the solar neighborhood as well as 
intrinsic (i.e. not caused by contamination eff'ects due to dynam- 
ics), and thus the one to be reproduced, this presents a serious 
problem for theoretical predictions. As already mentioned in the 
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Fig. 13. Fraction of all Fe produced by SNe la with the combined SD (with ci = 1) + DD model, with all double C-O WD mergers 
resulting in SNe la and a convective core mass increase of 10%, for a constant SFR (black) and for the two-infall model (gray). 



context of the G-dwarf distribution, almost no theoretical chem- 
ical evolution models are able to produce an AMR that is flat 
throughout most of the Galaxy's history, although again the ra- 
dial migration assumption may help i n this respect. H owever, an 
alternative AMR was derived by Ram irez et al.l (|2007|) under dif- 
ferent conditions on the star selection and especially the age un- 
certainty, as well as different stellar parameter zero points (espe- 
cially temperature, leading to different ages). This observational 
AMR does indeed show a rise in function of time, as is predicted 
by theoretical evolution models. However, it is likely that this 
sample consisting of stars suitable for high-resolution high S/N 
spectral analysis is not representative for the solar neighborhood, 
as it is based on a collection of spectroscopic data for a very 
limited number of stars, selected in a non-random fashion, but 
poorly constrained nonetheless (1. Ramirez, personal communi- 
cation, 2012). Nev ertheless, this AMR is a lso shown in Fig. |9l 
Wheth er or not the iHolmberg et al.l (l2007l) or ICasagrande et al.l 
(l20Tll) AMR are indeed representative for the solar neighbor- 
hood additionally depends on the validity of the used stellar 
isochrones. For completeness, it should also be noted that if 
one were to constru ct a G-dwarf metallic ity distribution from 
the stars selected by iRamirez et al.l (l2007h for the construction 
of their AMR, this G-dwarf metallicity d istribution would beex- 
tremely more spread out than the one bv Holmberg et al] (l2007h 
or ICasagrande et al. (2011). While such distribution would be 
questionable due to the mentioned non-representative nature of 
the dataset, it is of interest to note that to retrieve such a spreaded 
G-dwarf metallicity distribution in our predictions, the intrinsic 
scatter added to the time-dependent [Fe/H] must be taken very 
large, to the extent where it would correspond to the entire scat- 
ter observed in the AMR. 

The evolution of the SN la rate as a function of time is shown 
in Fig. [To] with the currently observed rate indicated. It shows an 
acceptable level of agreement with observation for either model, 
although it should be remembered from Fig. |7] that only for the 
flat SFR a good agreement in morphological shape was found 
for the G-dwarf metallicity distribution. 

Figures [TT] and [12] show the C and O abundance vari- 
ation obtained with the same model, again for both flat 
SFR and the two-infall mode l . Obs e rved C abu n dance s 
are taken from IChiappini et"an (l2003h . iReddv et alJ (l2003h . 



'Bensbv & Feltzing' (2006*) and 'Fabbian et al.' (2009'), O abun- 
danc es from Gratt on et al. (2003), Reddy et al. (2003) and 
Cal ura et al.1 (l2010h . While early [C/Fe]-values seem to lie 
slightly too high, and the decreasing trend between [Fe/H]=- 
0.5 and suggested by observations is not reproduced, the mod- 
els are certainly not incompatible with them. As mentioned be- 
fore, the early evolution ([Fe/H]<-1.5) of both elements can be 
greatly modified by making different assumptions about direct 
black hole formation and hypernova yields, but this does not af- 
fect the G-dwarf metallicity distribution which is the object of 
this study. 

As a final step, the fraction of Fe in the solar neighborhood 
originating from SNe la is investigated. Fig. [13] shows its evo- 
lution as a function of time, for the constant SFR, respectively 
two-infall model calculations. Between 65 and 70% of the Fe 
present in the solar neighborhood today was created by SNe la. 

3.3. Other considerations 

As mentioned in the Introductio n, concurrent with t his study 
a comparison has been made by iToonen et alJ (1201 3h between 
different binary population synthesis models. This exercise has 
brought to light a number of ingredients, other than intended pa- 
rameterizations, which may have a profound influence on the 
outcome of the evolution of a certain system. We will now inves- 
tigate their consequences for the results obtained so far. A first 
point is the wind mass loss during the AGB phase . Its yields as 
computed bv^ van den Hoek & GroenewegenldT99 7') are included 
in the code. However, we do not follow the angular momen- 
tum loss as a result of these winds in detail. Instead, the mass 
decrease from TAMS to post-mass transfer phase (i.e. the C-O 
core) and the resulting orbital separation variation is calculated 
in one step, thus in fact acting as if this mass is lost in the CE. Of 
course, this results in a different orbital separation for individual 
systems as in those codes which follow the AGB phase in detail 
throughout at the expense of more computation time. However, 
the effect of this simplification on the eventual DTDs (and thus 
also the G-dwarf metallicity distributions) has been found to be 
negligible. The same is true for the simplification of not follow- 
ing in detail the convective core mass evolution between systems 
filling their Roche lobe early and late on the AGB. 
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It is important to remark that the considered mass range of 
the IMF is of paramount importance to the obtained absolute rate 
of SNe la in theoretical studies. If the lower limit of the IMF is 
taken at about 0. 1 M© (as most IMFs actually prescribe), then a 
large portion of the total mass of a starburst is "locked" in low 
mass systems that do not have a chance to ever result in a SN la 
(or even return their mass to the interstellar medium in any way). 
If, however, one takes the lower limit at 1.0 M©, the contribution 
to the total mass of the starburst of those systems that are too 
light to ever be able to produce a SN la is significantly reduced, 
and the SN la rate expressed in SNuM (events per formed stellar 
mass) increases. It is well known that most theoretical studies 
produce absolute SN la rates that are up to an order of magnitude 
lower than the rates that are being observed. One could thus be 
tempted to use an adapted IMF range to alleviate this problem. 
However, this should not be done unless there are indications 
that the observationally studied populations have a low (or even 
zero) contribution of low mass stars, which is not the case. 



4. Conclusions 

We have tested the ability of the two most popular type la su- 
pernova progenitor scenarios, the single and double degener- 
ate channel, to theoretically reproduce the observed indicators 
of the chemical history of the solar neighborhood. In particular, 
the metallicity distribution of G-type dwarfs (i.e. the prevalence 
of such stars with respective iron content) is such an indicator, 
which is additionally critically affected by the type la supernova 
rate throughout the history of the Galaxy, and thus by the as- 
sumptions about their progenitors. This was done by the combi- 
nation of a population synthesis code including detailed binary 
evolution and a galactic chemical evolution model. We applied 
various assumptions and parameters concerning both binary and 
galactic evolution. Additionally, we updated and extended the 
study by Mennekens et al. (2010) concerning the ability of both 
progenitor scenarios to reproduce the observed type la super- 
nova delay time distribution (the evolution of the number of such 
events as a function of time after starburst) in passively evolving 
elliptical galaxies. 

We conclude that both the delay time distribution and the 
G-dwarf metallicity distribution point towards a significant con- 
tribution by both the single degenerate and double degenerate 
channel. The ability of either traditional progenitor scenario to 
be solely responsible for all type la supernova induced iron en- 
richment is ruled out. The best match with the observed G-dwarf 
metallicity distribution (and other galactic observables) is found 
when assuming that all double C-O white dwarf mergers explode 
as type la supernovae, as well as using a slightly larger convec- 
tive core mass to account for rotational effects. In addition, the 
critical dependence of both distributions on certain binary and 
galactic evolutionary processes, exactly those processes which 
are modeled with still uncertain parameters in population syn- 
thesis studies, is a way to learn more about these processes and 
thus further constrain these parameters. In particular, we find that 
most double degenerate type la supernova progenitors need to go 
through a quasi-conservative, stable Roche lobe overflow phase, 
followed by a common envelope evolution. On a galactic level, 
a satisfying reproduction of the G-dwarf metallicity distribution 
requires a high binary frequency and a star formation rate that is 
relatively constant in time. 

Acknowledgements. We thank Ivan Ramirez for a useful discussion regarding 
our comparison witli observational AMRs. 
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